Impact of Secondary non-Gaussianities on the Search for Primordial Non-Gaussianity 

with CMB Maps 
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When constraining the primordial non-Gaussianity parameter /nl with cosmic microwave back- 
ground anisotropy maps, the bias resulting from the covariance between primordial non-Gaussianity 
and secondary non-Gaussianities to the estimator of /nl is generally assumed to be negligible. We 
show that this assumption may not hold when attempting to measure the primordial non-Gaussianity 
out to angular scales below a few tens arcminutes with an experiment like Planck, especially if the 
primordial non-Gaussianity parameter is around the minimum detectability level with /nl between 
5 and 10. In future, it will be necessary to jointly estimate the combined primordial and secondary 
contributions to the CMB bispectrum and establish /nl by properly accounting for the confusion 
from secondary non-Gaussianities. 
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Introduction — The search for primordial non- 
Gaussianity with constraints on the non-Gaussianity pa- 
rameter /nl using cosmic microwave background (CMB) 
anisotropy maps is now an active topic in cosmol- 
ogy today [J [3, [1]. The 3-year Wilkinson Microwave 
Anisotropy Probe has allowed the constraint that — 54 < 
/nl < 114 at the 95% confidence level [l], though a more 
recent study claims a non-zero detection of primordial 
non-Gaussianity at the same 95% confidence level with 
26.9 < /nl < 146.7 This result, if correct, has 
significant cosmological implications since the expected 
value under standard inflationary models is /nl ^ 1 
B S B B E El, El , though alternative models of infla- 



tion, such as the ekpyrotic cosmology [ij, ll5| , generally 
predict a large primordial non-Gaussianity with /nl at 
few tens. 

Most studies that constrain /nl with CMB anisotropy 
maps make use of an estimator for /nl of the form jl6l , 
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when Bf^i^f,^ is the primordial bispectrum with the as- 
sumption that /nl = 1- Here, 5*'™ is a linear correction 
to account for issues related to the maps (such as the 
mask) and N is an overall normalization factor [3] ■ For 
the present discussion motivated from an analytical cal- 
culation, we can ignore the correction associated with 
5''" which involves imperfections in the data, such as 
due to the mask. We also assume all-sky data here. In 
equation ([2]), is the noise variance to the bispectrum 



In general 



U ^2 



A rSZ-k, 4 oISW- 



of the non-Gaussianity with an overall normalization 
given by 6ps and Bf^~i^ and B]^^i~'^ are additional fore- 
ground, secondary non-Gaussianities from the SZ and 
ISW effects correlating with CMB lensing These 
are certainly not all the non-Gaussian contributions to a 
CMB map. There are non-Gaussianities from ISW [2H, 
kinetic SZ/Ostriker-Vishniac [ll, and the SZ effect itself 
[H. We ignore ISW and kinetic SZ/OV related bispec- 
tra as they are small compared to SZ generated non- 
Gaussianities. The SZ-SZ-SZ bispectrum is significant at 
arcminute angular scales, but given the power-law shot- 
noise behavior of the SZ bispectrum when I < 1500, the 
SZ contribution to the bispectrum can be thought of as 
an additional correction to 6ps. The shot-noise behavior 
of the SZ effect is especially applicable for the SZ contri- 
bution during reionization associated with hot electrons 
in supernovae bubbles Compton-cooling off of the CMB 
2J|. Thus, we do not separately include the total SZ 
bispectrum as a separate non-Gaussianity here. 

When estimating /nl, it is usually assumed that 
°ihh ~ /wL^f"™ when estimating the primordial bis- 
pectrum. This allows an estimator for /nl through 
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with /nl — Spi-im/N, where the normalization N is sim- 
ply the summed term. The above assumption that only 
the primordial non-Gaussianity can be considered is gen- 
erally motivated by the fact that the covariance term 



associated with the mode overlap between B 
additional secondary contributions to B' 
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where Bf^, , 



is the shape 



when Ai = (6ps, ^sz, ^isw, ■••) is expected to be smaller 
than the dominant term from equation ([3]) [l[. Never- 
theless, an estimate of /nl only from equation ([3]) leads 
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FIG. 1: Absolute values of signal-to-noise ratio squared for 
the detection of primordial bispectrum (black lines) assuming 
/nl = 1 as a function of I3. The signal-to-noise ratio squared 
for WMAP and Planck are shown with dashed and dot-dashed 
line, respectively. The red, blue, and green lines show the 
confusion resulting from the covariance between primary and 
point source, primary and SZ-lensing, and primary and ISW- 
lensing bispectra, respectively. 



to a biased estimate because of the contributions from 
secondary anisotropics through equation 

While the CMB map contains a large number of sec- 
ondary non-Gaussian signals, in terms of the covariance 
related to the /nl measurement, what is necessary is 
not to account for all of these non-Gaussianities, but 
to account for non-Gaussianities with bispectrum shapes 
Bi-^i^i^ in {11,12,13) moment space that align with the 
shape of the primary bispectrum. In this respect, pre- 
vious calculations have suggested that the point-source 
bispectrum may be ignored [l| , but the ISW-lensing bis- 
pectrum must be accounted for the Planck analysis [25| . 

Including the SZ-lensing bispectrum, we find that 
while the assumption that the covariance from secondary 
anisotropics can be mostly ignored for an experiment like 
WMAP, it may be necessary to account for certain covari- 
ances when estimating /nl from a high resolution exper- 
iment like Planck, especially if the underlying primordial 
non-Gaussianity has a value around /nl between 5 and 10 
consistent with the minimum amplitude detectable with 
Planck. At the minimum detectability level of WMAP 
with /nl ~ 20, the secondary anisotropics involving both 
residual points sources and lensing correlations will bias 
/nl by a factor between 1.2 and 1.5 if primordial non- 
Gaussianity estimate is performed out to angular scales 
corresponding to ^ > 700. 

To reach these conclusions, we first calculated iJf'VT 
following Ref. [1| with the full radiation transfer function 
using a modified code of CMBFAST ^ for the standard 
fiat ACDM cosmological model consistent with WMAP 
with fib = 0.042, = 0.238, h = 0.732, n = 0.958, and 



T = 0.089. We verified our calculations are consistent 
with prior calculations in the literature. In Fig. 1, we 
show the the absolute value of the signal-to-noise square 
ratio for the primordial bispectrum (thick lines) and for 
the covariances between primary and secondary bispec- 
tra. The plotted quantity here involving d{S/N)^/d\nl3 
resembles the estimator S above, expect for the sum 
over I3 while keeping the sign (ignoring the sign changes 
lead to a higher bias as described in Ref. [25|). While 
the primordial calculation involves ^(i?P''™)^/cr^, the 
"signal-to-noise" square of the covariance follows from 
6ps X] BP"'^BP^/a'^, for example for the point-source con- 
fusion, and these confusions should not be interpreted 
simply as the signal-to-noise ratio square to detect any of 
these secondary bispectra directly from the GMB maps. 

Instead of squared signal-to-noise ratios, to highlight 
the bias introduced to /nl when the estimator ignores 
secondary non-Gaussianity covariances, we calculated 
/nl = /nl + /bias where /bias is the bias that is gen- 
erated artificially by the correlation of modes between 
the primordial bispectrum and secondary bispectra. To 
properly normalize the relative contribution from sec- 
ondary non-Gaussianities, we assume normalizations for 
the point-source bispectrum consistent with WMAP with 
&ps = 3 X 10~^^, consistent with Q+V+W residual fore- 
ground [4], and Planck with = 5 x 10"^''. The 
&ps value for Planck is slightly higher than the values 
routinely quoted in the literature for unresolved radio 
sources in Planck high resolution maps, but this is due 
to the fact that we believe bps includes additional con- 
tributions such as from the SZ-SZ-SZ bispectrum from 
both clusters at low redshifts and supernovae halos dur- 
ing reionization with a power-law shot-noise spectrum 
when / < 1500. For the ISW-lensing and SZ-lensing bis- 
pectrum, we follow the calculation of Ref. and gen- 
erate the SZ contribution and the SZ correlation with 
dark matter halos responsible for lensing of the CMB 
using the halo model 23]. To account for an overall un- 
certainty and the variation in SZ and ISW amplitudes 
we have introduced an overall amplitude A^z and Ajsw 
respectively. Finally, to illustrate our results, we assume 
/nl consistent with roughly the minimum detectable pri- 
mordial non-Gaussianity with WMAP and Planck with 
/nl = 20 and 5, respectively. As we find later, the dom- 
inant confusion is from lensing bispectra and not from 
point sources. 

We summarize our results in Fig. 2, where we plot 
/n'l which can be thought of as the total primordial 
non-Gaussianity parameter that one will extract with 
the above estimator for /nl when no attempt has been 
made to separate out the confusion from secondary 
anisotropics. For the most part, the bias is negligible 
and becomes only important when I > 500. For WMAP, 
shown with a dashed line in Fig. 2 with the assumption 
that /nl = 20 if non-Gaussianity measurements are at- 
tempted out to / > 700, capturing basically all informa- 
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FIG. 2: The maximum non-Gaussianity measured with an 
optimal estimator for the primordial bispectrum /nl , which 
includes the true underlying primordial non-Gaussianity with 
/nl as labeled on the figure and the bias correction coming 
from the unaccounted secondary anisotropies. The bias is 
generally small and non existing if primordial non-Gaussianity 
measurements are limited out to I < 500, but depending on 
the value of /nl and the residual point source contamination, 
the correction is generally a factor of 1.5 to 2. If /nl ^ 
10, for Planck, it is necessary to account for secondary non- 
Gaussianities properly. 

tion in WMAP maps, then one finds a bias between a fac- 
tor of 1.5 to 2 if /nl ~ 20. If /nl > 30, then the relative 
contribution from secondary non-Gaussianities are sub- 
dominant compared to the primordial non-Gaussianity. 
Alternatively, if WMAP data are used to constrain that 
/nl < 30, then such a constraint must account for the 
covariances from secondary non-Gaussianities, especially 
those involving CMB lensing. 

With Planck, non-Gaussianity estimates can be ex- 
tended to ^max ~ 2000, but at such small angular scales, 
one finds a bias higher by a factor of more than 2 rela- 
tive to the lowest value of /nl that can be reached with 
Planck (dot-dashed line). In return, if Planck data were 
to constrain /nl to be below ~ 20, then such a con- 
straint must account for the confusion from secondary 
anisotropies to the "optimal" estimator of /nl, since lens- 
ing non-Gaussianities produce a correction to /nl with 
/bias ~ 10. 

To account for secondary non-Gaussianities, one can 
modify existing "optimal estimators" for /nl and jointly 
fit for both the primordial non-Gaussianity and the sec- 
ondary non-Gaussianities through a series of estimators 
Sa where a denotes the non-Gaussianity of interest with 

So. = N^jjKfj , (5) 

where 

111 {h,h,h) 



and refers to the set of non-Gaussianity parameters: 
(/nl, bps, ^sz, ^isw)- This method assumes that one has 
a good model for (Zi, Z2, ^3) dependence of secondary bis- 
pectra -B/^/^jg- Even if the point source covariance is 
small, the amplitude of the point source confusion is gen- 
erally unknown. Moreover, at / < 1500 many secondary 
bispectra such as the SZ effect has a power-law behav- 
ior similar to the bispectrum of point sources. Thus, it 
would be necessary to determine the amplitude 6ps from 
a joint fit. 

Our suggestion that an accounting of secondary 
anisotropies is necessary for primordial non-Gaussian 
measurement is different from the general assumption in 
the literature that one can simply ignore the covariance 
between primordial and secondary non-Gaussianities. 
This partly comes through, for example, the suggestion 
that primordial and point-source bispectra are orthogo- 
nal following results from an exercise that involved jointly 
measuring non-Gaussian amplitudes /nl and 6ps using a 
set of simulated maps in Ref. Q to study if there are 
biases in the estimators. However, this study used simu- 
lated non-Gaussian maps that did not include any point 
sources with bps = 0. This sets the covariance to be zero 
and we believe this may have led to the wrong conclu- 
sion that there is no bias in the optimal estimator for 
/nl from unresolved point sources, though such a bias is 
expected to be small, but non-negligible if /nl ^ 1- Our 
conclusions are consistent with some of the observations 
in Ref. [H. 

Here we have considered the confusion from secondary 
non-Gaussianities such as point sources and those gen- 
erated by CMB lensing. Additional contributions to the 
bispectrum exist with correlations between SZ, ISW and 
point sources as they all trace the same large-scale struc- 
ture at low redshifts. Previous studies using the halo 
model to describe the non-linear density field have shown 
correlations such as between SZ and radio sources to be 
small f23( , but since the bispectra in these cases are of the 
form SZ-PS-PS"', these bispectra may have a multipolar 
dependence in (^1,^2,^3) that is more aligned with the 
CMB primary bispectrum. In an upcoming paper (27l | 
we will discuss the impact of such foreground bispectra 
due to correlations between CMB secondary anisotropies 
and point sources. While our discussion has concentrated 
on a momentum independent non-Gaussianity parameter 
/nl, or the so-called local type associated squeezed tri- 
angles, it is easy to generalize the calculation for more 
complex descriptions of /NL(ki, k2, ks) [2§|. Due to dif- 
ferences in mode overlap, the exact momentum depen- 
dence will change the covariance contributions and the 
impact of secondary non-Gaussianities will be different 
between attempts to measure local /nl and, for exam- 
ple, equilateral /nl- 

Based on our calculations on the covariance between 
lensing and primary bispectra we have suggested a po- 
tential confusion for /nl measurement in Planck data. It 
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is unlikely that our observation on the importance of sec- 
ondary non-Gaussianities changes any of the current con- 
straints on the non-Gaussianity parameter with WMAP 
data given that they mostly lead to /nl ^ 100 roughly. 
The secondary non-Gaussianities, however, could im- 
pact the significance of any detections of primordial non- 
Gaussianity, especially if the detection is marginally dif- 
ferent from zero [1]. For such studies, the exact signif- 
icant of the detection should include an accounting of 
the secondary non-Gaussianity and the overlap with pri- 
mordial bispectrum in the "optimal" estimator used to 
establish /nl- 

We thank Eiichiro Koniatsu for a helpful communica- 
tion. This work was supported by NSF CAREER AST- 
0645427. We acknowledge the use of CMBFAST by Uros 
Seljak and Matias Zaldarriaga [1^. 



[1] 
[2] 

[3] 
[4] 
[5] 
[6] 

[7] 



[9] 
[10] 

[11] 



E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 
(2001) ar Xiv:astro- ph/0005036 . 

M. Liguori, F. K. Hansen, E. Komatsu, S. Matar- 
rese and A. Riotto, Phys. Rev. D 73, 043505 (2006) 
arXiv:astro-ph/0509098 . 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. 

Suppl. 148, 119 (2003) arXiv:astro-ph/0302223 . 

D. N. Spergel et al, Astroph^JTSuppL 170, 377 (2007) 

|arXiv:astro-ph/0603449, 

A. P. S. Yadav and B. D. Wandelt, [arXiv:astro- 
ph0712.1148]. 

D. S. Salopek and J. R. Bond, Phys. Rev. D 42, 3936 
(1990); ibid. 43, 1005 (1991) 

T. Falk, R. Rangarajan and M. Srednicki, Astrophys. J. 
403, LI (1993) 

A. Gangui, F. Lucchin, S. Matarrese and S. MoUerach, 
Astrophys. J. 430, 447 (1994) 

T. Pyne and S. M. Carroll, Phys. Rev. D 53, 2920 (1996) 

J. M. Maldacena, JHEP 0305, 013 (2003). 

V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, 



[12; 

[is: 

[14 

[is: 

[16 
[17 

[is: 

[19 
[20 
[21 

[22: 

[23: 

[24: 

[25: 

[26: 

[27 

[28: 



Nucl. Phys. B667, 119 (2003), [arXiv:astroph-/0209156] 
N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. Lett. 
93, 231301 (2004) [arXiv:astro-ph/0407505| . 
N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, 
Phys. Rept. 402 (2004) 103-266 

E. L Buchbinder, ,L Khoury and B. A. Ovrut, 
larXiv:0710.5172 [hep-th]. 

J. L. Lehners and P. J. Steinhardt. larXiv:0712. 37791 [hep- 
th]. 

E. Komatsu, D. N. Spergel and B. D. Wandelt, Astro- 
phys. J. 634, 14 (2005) arXiv:astro-ph/0305189' . 

P. Creminelli, A. Nicolis, L Senatore, M. Tegmark and 
M. Zaldarriaga, Journal of Cosmology and Astro-Particle 
Physics 5, 4 (2006), arXiv:astro-ph/0509029 
A. P. S. Yadav, E. Komatsu, B. D. Wandelt, M. Liguori, 

F. K. Hansen and S. Matarrese, ,arXiv:0711.4933l [astro- 
ph]. 

A. R. Cooray and W. Hu, Astrophys. J. 534, 533 (2000) 
arXiv;astro-ph/9910397 . 

D. M. Goldberg and'D. N . Spergel, Phys. Rev. 
103002 (1999) arXiv:astro-ph/9811251,. 
A. 



D 59, 



Cooray, Phys. Rev. 

jarXiv:astro-ph/0109162' ; 
D. N. Spergel, Phys. Rev, 
[arXiv:astro-ph/0 108179 . 
A. Cooray, Phys. Rev. 



D 65, 



D 65, 



083518 
Verde 
043007 



(2002) 
and 
(2002) 



D 



64, 063514 (2001) 
[arXiv: astro-ph/0105063 ; P. G. Castro, Phys. Rev. 
D"67, 044039 (2004) [Erratum- ibid. D 70, 049902 
(2004)] arXiv:astro-ph/0212500 . 

A. Cooray and R. K. Sheth, Phys. Rept. 372, 1 (2002) 
arXiv:astro-ph/0206508 ; A. Cooray^ P hys. Rev. D 62, 
103506 (2000) arXiv:astro -ph76605287| . 
S. P. Oh, A. Cooray and M. Kamionkowski, 
Mon. Not. Roy. A stron. Soc. 342, L20 (2003) 
[arXiv:astro- ph/0303007 
K. M. Smith 
arXiv:astro-ph/0612571 
U. Seljak and M. Zaldarriaga, 
(1996) arXiv:as tro-ph/9603033 
D. Sarkar et al. in preparation. 
D. Babich, P. Creminelli and 



and 



M. 



Zaldarriaga, 



Astrophys. J. 469, 437 



0408, 009 (2004) [arXiv:astro-ph/0405356] . 



M. Zaldarriaga, JCAP 



